Primordial gravity waves fossils and their use in testing inflation 
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A new efTect is described by which primordial gravity waves leave a permanent signature in the 
large scale structure of the Universe. The efTect occurs at second order in perturbation theory and 
is sensitive to the order in which perturbations on different scales are generated. We derive general 
forecasts for the detectability of the effect with future experiments, and consider observations of the 
pre-reionization gas through the 21cm line. It is found that the Square Kilometre Array will not 
be competitive with current cosmic microwave background constraints on primordial gravity waves 
from inflation. However, a more futuristic experiment could, through this effect, provide the highest 
ultimate sensitivity to tensor modes and possibly even measure the tensor spectral index. It is thus 
a potentially quantitative probe of the inflationary paradigm. 



Introduction. — It has been proposed that redshifted 
21cm radiation, from the spin flip transition in neutral 
hydrogen, might be a powerful probe of the early uni- 
verse. The era before the first luminous objects reionized 
the universe-around redshift 10-contains most of the ob- 
servable volume of the universe, and 21 cm radiation is 
the only known probe of these so called dark ages (see 
Furlanetto et al. fl] for a review). The density of the 
hydrogen could be mapped in 3D analogous to how the 
cosmic microwave background (CMB) is mapped in 2D. 
The wealth of obtainable statistical information may al- 
low for the detection of many subtle effects which could 
probe the early universe. In particular, the primordial 
gravity wave background, also referred to as tensor per- 
turbations, are of considerable cosmological interest. 

Inflation robustly predicts the production of tensor 
perturbations with a nearly scale-free spectrum, however, 
their amplitude is essentially unconstrained theoretically. 
The amplitude of the tensor power spectrum is quanti- 
fied by r, the tensor to scalar ratio. The current upper 
limit is r < 0.24 at 95% confidence however upcom- 
ing CMB measurements will be sensitive down to r of a 
few percent The current limits on r correspond to 
characteristic primordial shear on the order of 10^"' per 
logarithmic interval of wavenumber. 

Several probes of gravity waves using the pre- 
reionization 21 cm signal have been proposed. These in- 
clude polarization Q and redshift space distortions 
Dodelson et al. @ considered the weak lensing signature 
of gravity waves and found that the signal is sensitive to 
the so called metric shear. This is closely related to the 
present work. 

Here we describe a mechanism by which primordial 
gravitational waves may leave an imprint in the statistics 
of the large scale structure (LSS) of the universe. This 
signature becomes observable when the gravity wave en- 
ters the horizon and begins to decay. 

Mechanism. — In the following, Greek indices run from 
to 3 and lower case Latins from 1 to 3. Latin indices are 
always raised and lowered with Kronecker deltas. Com- 



mas denote partial derivatives, and an over-dot (#) rep- 
resents a derivative with respect to the cosmological con- 
formal time. Finally, we adopt a mostly positive metric 
signature (—1, 1, 1, 1). 

We start with an inflating universe with some distribu- 
tion of previously generated tensor modes that are now 
super horizon. Scalar, vector and smaller scale tensor 
modes may exist but their contribution to the metric is 
ignored. The line element is given by 



(1) 



where a is the scale factor, rj the conformal time and 
a spatially flat background geometry has been assumed. 
The metric perturbations hij are assumed to be trans- 
verse and traceless and thus contain only tensor modes. 
The elements of hij are also assumed to be small such 
that only leading order terms need be retained. The as- 
sumption that all tensor modes under consideration are 
super horizon implies that kh ^ d/a, where kh denotes 
the wave numbers of tensor modes. The frame in which 
the line element takes the form in Eq. [1] will hereafter be 
referred to as the cosmological frame (CF). 

By the equivalence principle, it is possible to perform 
a coordinate transformation such that the space-time ap- 
pears locally Minkowski at a point. New coordinates are 
deflned in which the tensor modes are gauged away at 
the origin: 



(2) 



where the elements hoa are taken to be zero. The metric 
now takes the form (up to flrst order in hij) 

ds^ = [—di]^ + Sijdx^dx-' — x"^ dah^cdx"' dx^'\ . (3) 

This frame will be loosely referred to as the locally Fried- 
mann frame (LFF), because in these coordinates the met- 
ric is locally that of an unperturbed FLRW Universe. We 
will give quantities in these coordinates a tilde (#) to 
distinguish them from their counterparts in the CF. It 
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is seen from Eq. |3] that the local effects of gravity waves 
are suppressed not only by the smallness of hij but also 
by kfi/k where k = L^^ and L is some length scale of 
interest. This will be important in justifying some later 
assumptions. Note that for super horizon gravity waves, 
temporal derivatives are much smaller than spacial ones. 

On small scales, inflation generates scalar perturba- 
tions which are then carried to larger scales by the expan- 
sion. By the equivalence principle, physical processes on 
small scales can not know about the long wavelength ten- 
sor modes. As such these small scale scalar modes must 
be uncorrelated with the long wavelength tensor modes. 
We assume statistical homogeneity and isotropy in the 
LFF as would be expected from inflation. The power 
spectrum of scalar perturbations can then be written as 
a function of only the magnitude of the wave number, 
i.e., P{ka) — P{k). This applies only within the local 
patch near the point where the tensor mode was gaged 
away. The average in the definition of the scalar power 
spectrum is over realizations of the scalar map, but not 
the tensor map. 

In the CF, the isotropy is broken. Transforming back 
to cosmological coordinates maps ki — > ki~kjh^-' /2. The 
power spectrum becomes sheared: 

P{ka) = P{k) - + Oi'^h,,) + 0{h.,,\ (4) 

If the metric perturbations are not assumed to be trace- 
less, the right hand side of this equation gains an ad- 
ditional term proportional to this trace. This deviation 
from isotropy is not observable since any possible obser- 
vation would take place in the LFF. 

It is noted that the leading order correction to CF 
power spectrum is not suppressed by kh/k. It is therefore 
not expected that the residual terms in the LFF metric 
(Eq. ISj can break isotropy to undo CF anisotropy. How- 
ever it was the CF in which the power spectrum should 
be isotropic, then there would be observable anisotropy 
in the LFF. This would be a violation of the equivalence 
principle, since an experiment local in both space and 
time would be able to detect the super horizon tensor 
modes by measuring the power spectrum of the locally 
generated scalar perturbations. 

We would now like to evolve the system to some later 
time when observations can be made. Ignoring the in- 
ternal dynamics of the scalar perturbations, we solve for 
their evolution as if they were embedded in a sea of test 
particles. This is trivial since an object at coordinate rest 
in the CF will remain at rest for any time dependence of 
hij (this is true at all orders). At some point well after 
inflation, when the universe is in its deceleration phase, 
the horizon will become larger than the length scale of 
the tensor modes. The tensor modes will then decay by 
redshifting, and after some period of time the metric per- 
turbations hij become negligible. The CF and LFF then 
become equivalent and both correspond to the frame in 



which observations can be made. The distribution of test 
particles is the same as it initially was in the CF. As such, 
the initially physically isotropic power spectrum now con- 
tains a measurable local anisotropy given by Eq. |4l The 
values of the initial metric perturbations can be deter- 
mined by measuring this distortion at any time in the 
future, constituting a fossil of the initial tensor modes. 

The scalar perturbations remain Gaussian but be- 
come non-stationary, and the trispectrum gains the cor- 
responding terms. This is analogous to the apparent dis- 
tortions expected in the CMB and 21 cm fields induced 
by gravitational lensing. Similarly the bispectra of mixed 
scalars and tensors were calculated in Maldacena Q , em- 
ploying similar methodology to that presented here. 

The effect described here is a second order perturba- 
tion theory effect, in that it is a small effect due to tensor 
modes on the already small scalar perturbations. This 
coupling occurs in the initial conditions, not between the 
dynamics of the scalars and tensors. The simple argu- 
ment presented above avoided the complication of a full 
second order calculation, but it is expected that such cal- 
culations would yield the same results. Specifically, an 
expression agreeing with Eq. Slto relevant order, was 
derived in Giddings and Sloth [8|, Eq. 4.5] as part of a 
longer calculation. 

Tests of inflation. — The above arguments relied on 
perturbations on large scales being generated before per- 
turbations on small scales. This is the case in any con- 
ceivable model of inflation, however it is not be the case 
in all scenarios. As an illustrative example, in the cos- 
mic defect scenario perturbations are generated on small 
scales and then causally transported to larger scales as 
the universe evolves. It is argued that in this scenario, 
tensor perturbations leave no fossils. A detection of pri- 
mordial tensors by another means (CMB B-modes for 
example) with an observed lack of the corresponding fos- 
sils would provide a serious challenge to inflation. 

The most speciflc prediction of single fleld inflation is 
the power spectrum of tensor modes, defined by 

{2^f5{ka - k'^)Ph{ka) = {h,j{kaW^k'^)). (5) 

Given the amplitude of the scalar power spectrum A^, 
the tensor power spectrum is fixed by a single parameter, 
the tensor to scalar ratio r. The shape of the spectrum 
is then nearly scale-free: 



We follow the WMAP conventions for defining Ph, As 
and r [§]. The spectral index fixed by the consistency 
relation, nt — —r/8 [3]. The pivot scale is taken to be 
ko = 0.002 Mpc"^ and we assume the WMAP7 central 
value for As of 2.46 x 10"^. 

Because r is likely small, any deviation from a scale- 
free spectrum will be difficult to measure, making the 
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verification of the consistency relation correspondingly 
difficult. The CMB is sensitive primarily to large scale 
tensor modes, with smaller scale modes having decayed 
by recombination. Cosmic variance and lensing contam- 
ination will likely prevent a measurement of rit from the 
CMB, unless the lensing can be cleaned from the sig- 
nal Conversely, the amplitude of the fossil signal 
does not decay as the universe expands. It may thus be 
possible to make a measurement of the spectral index, 
provided r is sufficiently large. 

Statistical detection in LSS. — In practice, the tensor 
gravity wave fossils could be reconstructed by applying 
quadratic estimators to the density field. Aside from the 
increased dimensionality, this is identical to the manner 
in which lensing shear is reconstructed 12, 1^. Rather 
than considering the statistics of such estimators, here 
we follow a simpler line of reasoning to approximate the 
accuracy to which the tensor parameter can be measured. 

We begin by asking how well a long wavelength, tensor 
mode can be reconstructed from its effects on the scalar 
power spectrum (Eq. |4]). The metric perturbations are 
assumed to be spatially constant and take the form 

h^J ^ h+e+{z) + h^ef^iz) (7) 

where and e]!- are the polarization tensors and the z 
direction of propagation is chosen for convenience. The 
uncertainty on the scalar power spectrum is 



[AP{ka)f =2[P{ka)+N\ 



{S) 



where N is the noise power. We use a Fisher Matrix 
analysis to sum this information over all ka to determine 
the corresponding uncertainty on the shear and hx- 
Assuming an experiment whose noise is sub dominant 
to sample variance (N P), the resulting variance is 
inversely proportional to the number of modes surveyed: 
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(9) 



where h stands for either /i+ or hy (the superscript C 
indicates that the formula applies for spatially constant 
h), V is the volume of the survey and k^ax is set by the 
resolution of the survey. The constant of proportionality 
depends on the shape of the unsheared power spectrum 
P(fc), but to within a few tens of percent it is unity. 
21 cm emission will be difficult to observe on large scales 
[l|, however it is small scales that dominate the number 
of modes and thus the reconstruction. It is only the co- 
herence of small scale anisotropy that must be measured 
on large scales. 

Given the reconstruction uncertainty on a spatially 
constant shear, and the fact that reconstruction noise is 
scale independent (white) 12] , the noise power spectrum 
for spatially varying tensor modes is then 



Nh = 4T/(A/i^)2 = 4 



27r 



(10) 



The factor of four comes from the definition of the power 
spectrum in Eq. [SJ noting that (hijh^^) = 4(/i^). 

We now sum over ka^ to determine the signal to noise 
as a function of tensor power spectrum amplitude r. The 
signal to noise ratio squared is then 



rkupp 



p2 



2{Nh + Phf 
'-^ dkk'^ Pjik) 

. 2^2 (iV, + p,)2 



(11) 
(12) 



It is seen from the redness of the spectrum P^ (Eq. |6]) 
that the result is completely independent of the upper 
limit of integration. The same redness makes the final 
result extremely sensitive to the lower limit. As described 
above, the fossil of a primordial tensor mode can only 
be observed once the mode has decayed. This begins 
to happen when the scale of the gravity wave becomes 
comparable to the horizon scale, and as such, the largest 
scale observable mode has wavelength kiower ~ o-H. 

For an initial detection, we assume that noise domi- 
nates sample variance at each ka, i.e., Nh ^ P^. Setting 
the signal to noise ratio to be 2, for a 95% confidence 
detection, yields a minimum detectable amplitude of 
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(13) 



where Vh = (aH)^^. 

While the observability of 21 cm radiation depends on 
the reionizaton model, one regime in which a strong 
signal may exist is near redshift 15 [1[. The planned 
Square Kilometer Array (SKA) will aim to probe this 
era with 10 km baselines [14| . Assuming a survey volume 
of 200 (Gpc//i)'^ and a noiseless measurement, the limit 
on r achievable with SKA will be 



7.3 



1.2Mpc//i 



^ '200 (Gpc//i)3 3.3 (Gpc//i)3T 



V 



Vh 



(14) 

While this constraint is not competitive with current con- 
straints from the CMB, it is a strong function of the 
resolution of the experiment. The Low Frequency Array 
(LOFAR) for instance, has baselines extending to 400 km. 
However LOFAR will not have the sensitivity to probe 



the dark ages [15|. It is the physical shear due to grav- 



ity waves at the source that is being measured, and all 
light propagation effects, such as the lensing considered 
in Dodelson et al. @, have been ignored. 



^ From this point forward, ka will refer to the wave number of a 
tensor mode, not a scalar mode. The exception will be kmax 
which is the smallest scale at which a scalar can be resolved. 
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FIG. 1: Primordial tensor power spectrum obeying the con- 
sistency relation for r = 0.1. The solid line is the tensor 
power spectrum. Error bars represent the reconstruction un- 
certainty on the binned power spectrum for a perfect exper- 
iment, surveying 200(Gpc//i)^ and resolving scalar modes 
down to kmax = 168/i/Mpc. The dashed, nearly vertical, 
line is the reconstruction noise power. The non-zero slope of 
the solid line is the deviation from scale-free. 



Similar arguments are used to find the achievable er- 
ror on the spectral index n^. Properly considering the 
degeneracy with r, the error on rit is: 



Ant = F 



2tt 



1 



-I 1/2 



rA,V 



(15) 



where F is a function of the combination of parameters 
VH/{kl^axrAs). In the limit that Ph{k = aH) > Nh, 
which is the limit in which a measurement of nt is 
possible, F is approximately 6. Assuming the same 
volume and redshift as above, and that r = 0.1, the 
consistency relation is tested at the 2 sigma level for 
kmax — 168/i/Mpc. The tensor power spectrum and er- 
ror bars for this scenario are shown in Fig. [1] 

Such a measurement is very futuristic indeed, requiring 
a nearly filled array with greater than thousand kilometre 
baselines. Note that such an experiment would be sensi- 
tive to r down to the 10~^ level. Also, higher redshifts 
contain even more information, though their observation 
is technically more challenging. 

Discussion. — Aside from the technical challenge of 
mapping the 21cm signal over hundreds of cubic giga- 
parsecs and down to scales smaller than a megaparsec, 
there may be other competing effects that could hinder 
a detection. Of primary concern is weak lensing which 
also shears observed structures, creating apparent local 
anisotropics. The weak lensing shear is of order a few 
percent, and is thus many orders of magnitude greater 
than gravity wave shear. However, the 3D map of grav- 
ity wave shear will be transverse, transforming intrinsi- 
cally as a tensor. To linear order, the lensing pattern is 
the gradient of a scalar. Even at higher order, lensing 
always maps one point in space to another and is thus at 
most vector like. This test does not exist for the CMB or 



lensing due to the lower dimensionality of these probes. 

Also of concern is the preservation of the anisotropy on 
small scales. The scale corresponding to fc = 168/i/Mpc 
is still larger than the Jeans length at these redshifts, and 
as such hydrogen should trace the dark matter. However, 
the evolution of scalar perturbations is mildly nonlin- 
ear, and it is possible that this evolution will erase the 
anisotropy. Detailed analysis of the nonlinear erasure of 
the anisotropy is deferred to future investigation. 

There has been much recent interest in searching for 
anisotropy, and this has some implications for the fossil 
signal. The constraints on quadrupolar isotropy in LSS 
by PuUen and Hirata [l6[ should already imply a weak 
constraint at the r < 10^ level. Constraints from the 
CMB are not relevant however, since modes spanning 
the surface of last scatter remain super horizon today. 

CMB B-modes will be the most sensitive probe of pri- 
mordial gravity waves in the next generation of experi- 
ments. However, fossils may eventually be sensitive well 
below the limits of the CMB. 
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